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ABSTRACT 

We  present  the  theory  and  multiphysics  simulation  of 
a  novel  type  of  photoelastic  modulator  (PEM)  for  use  in 
ultra-high  performance  FT  spectrometers  with  applica¬ 
tions  in  the  optical  recognition  of  explosives  and  CB  war¬ 
fare  agents.  This  new  technology  is  based  on  the  strain 
birefringence  induced  along  a  bar-shaped  optical  element 
by  a  phased  array  of  piezoelectric  transducers  (PZT).  Ap¬ 
propriate  phase  and  amplitude  relationships  among  the 
array  PZTs  ensure  that  induced  birefringence  oscillations 
add  constructively  and  result  in  high  total  retardation. 
When  used  in  PEM/FT  spectrometers  (Buican,  2006;  Bui¬ 
can  and  Carrieri,  2004,  2006),  monolithic  virtual- 
stack/phased-array  (VSPA)  PEMs  (Buican  2008)  will 
provide  ultra-high  scan  rates  (to  >105  interferograms/s); 
broad  spectral  range;  high  spectral  resolution;  and  high 
light  throughput. 

1.  INTRODUCTION 

For  historical  reasons,  commercial  off-the-shelf 
(COTS)  PEMs  have  been  designed  for  low  retardation 
amplitudes;  as  a  consequence,  higher  resolution  PEM/FT 
spectrometers  must  employ  stacked  PEMs  (up  to  ten  in 
one  design)  (Buican  and  Carrieri  2004).  However,  the 
very  large  number  of  optical  surfaces  in  such  a  system 
results  in  significant  light  losses.  We  cut  through  the  Gor¬ 
dian  knot  of  the  resolution-throughput  compromise  with  a 
PEM  devised  specifically  for  high-resolution,  high- 
throughput  PEM/FT  spectrometry.  Such  a  PEM  has  a 
single  bar-shaped  optical  element  in  lieu  of  a  stack  of 
conventional  PEMs,  with  multiple  PZTs  attached  to  the 
bar  at  appropriate  intervals.  As  light  travels  along  the  bar, 
local  retardations  add  up  to  total  values  comparable  to 
those  for  stacks  of  phase-locked  individual  PEMs.  The 
bar- shaped  PEM  can  thus  be  regarded  as  a  virtual  stack  of 
conventional  PEMs. 

The  difficulty  with  such  a  system  is  that  elastic  waves 
generated  by  individual  PZTs  propagate  not  only  across 
the  thickness  of  the  bar  but  also  along  its  length.  As  the 
associated  stress  and  consequent  birefringence  change 


sign  every  half- wavelength  of  the  longitudinally  propagat¬ 
ing  elastic  waves,  the  total  retardation  contributed  by  each 
PZT  will  oscillate  rather  than  monotonically  increase  with 
increasing  bar  length.  Furthermore,  overlapping  elastic 
waves  generated  by  neighboring  PZTs  have  differing 
phases,  thus  reducing  the  amplitude  of  the  associated  re¬ 
tardation.  These  factors  tend  to  defeat  the  very  purpose  of 
the  bar-shaped  optical  element  and  of  the  multiple  PZTs 
attached  to  it. 

We  present  in  this  paper  a  conceptual  framework  for 
solving  the  longitudinal  propagation  problem  and  under¬ 
standing  VSPA-PEM  operation.  This  framework  rests  on 
the  basic  idea  that  elastic  wave  propagation  along  the  axis 
of  the  bar  can  be  confined  through  the  use  of  “active  mir¬ 
rors”  consisting  of  additional  PZTs  driven  at  appropriate 
amplitudes  and  phases.  Furthermore,  we  show  that  a  suit¬ 
ably  driven  phased  array  of  PZTs  placed  along  the  PEM 
bar  behaves  like  a  stack  of  mechanically  isolated  conven¬ 
tional  PEMs — hence  the  virtual-stack/phased  array,  or 
VSPA,  name  of  the  technology.  Being  isolated,  these 
stacked  “virtual”  PEMs  can  be  driven  in  phase  and  thus 
constructively  contribute  to  a  total  retardation  which  is  the 
scalar  sum  of  the  individual  virtual  PEM  amplitudes. 

Finally,  we  discuss  the  strategic  advantage  offered  by 
VSPA-based  PEM/FT  as  the  enabling  technology  for  sev¬ 
eral  novel  techniques  for  (i)  the  detection  of  trace 
amounts  of  explosives  on  surfaces;  and  (ii)  the  early 
warning  of  dangerous  environmental  CB  threats  evolving 
on  the  battlefield  or  in  the  homeland. 

2.  PEM/FT  SPECTROMETRY  BACKGROUND 

PEM/FT  spectrometers  are  based  on  birefringence  in¬ 
terferometers,  which  introduce  a  variable  retardation  be¬ 
tween  interfering  waves  separated  by  their  mutually  or¬ 
thogonal  directions  of  polarization.  The  retardation  is 
proportional  to  the  instantaneous  birefringence  of  the 
modulator  and  can  be  varied  in  time  in  order  to  scan  an 
interferogram.  PEM-based  FT  spectrometers  cover  broad 
spectral  ranges  (UV  to  NIR  with  fused  silica  optics,  and 
NIR  to  mid-IR  with  ZnSe  optics).  In  its  simplest  imple- 
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mentation  (Fig.  la),  a  birefringence  interferometer  con¬ 
sists  of  a  birefringent  modulator  placed  between  two  pola¬ 
rizers  (Buican  and  Martin,  1990).  The  optical  axes  of  the 
birefringent  modulator  are  at  45°  relative  to  the  direction 
of  the  polarizers. 


Fig.  1.  PEM/FT  spectrometry,  (a)  The  basic  birefringence  interferome¬ 
ter.  Ejn,  Eout — input  and  output  electric  vectors;  M — birefringence  modu¬ 
lator;  nx,  ny — principal  indices  of  refraction  of  the  birefringence  modula¬ 
tor  and  the  corresponding  optical  axes;  and  Pb  P2 — polarizers  and  their 
optical  axes.  If  M  is  a  single  resonator  PEM,  then  this  is  also  a  schematic 
illustration  of  a  first-generation  PEM/FT  spectrometer,  (b)  Second- 
generation  PEM/FTIR  spectrometer  prototype  design  based  on  actively 
tuned,  stacked  conventional  PEMs:  A — entrance  aperture;  D — detector; 
B — beam  splitters;  M — mirrors;  CD — calibration  detector;  P — 
polarizers;  CS — calibration  source;  PEMs — stacked  photoelastic  mod¬ 
ulators.  (c)  Conceptual  design  of  a  third-generation  PEM/FTIR  spectro¬ 
meter  based  on  a  bar-shaped  virtual-stack  PEM,  driven  by  a  phased 
array  of  PZTs.  Components  are  as  labeled  in  (b)  and  OE — optical  ele¬ 
ment. 

First  generation  PEM/FT  spectrometers  use  single 
PEMs  (or  multiple  PEMs  without  special  driving  and  con¬ 
trol  systems  characteristic  of  the  second  generation).  Typ¬ 
ically,  a  first  generation  PEM/FT  spectrometer  is  mainly 
suited  for  low  spectral  resolution  applications  such  as  the 
analysis  of  molecular  fluorescence  (Buican,  1985). 

Second  generation  PEM/FT  spectrometers  (Fig.  lb) 
address  this  limitation  by  (i)  increasing  PEM  driving  am¬ 
plitude  and  (ii)  stacking  multiple  PEMs.  Both  aspects  of 
the  technology  require  dedicated  driving  and  control  elec¬ 
tronics  that  make  the  stacked  PEMs  operate  in  phase  by 
actively  and  continuously  tuning  them  to  a  common  reso¬ 
nant  frequency  (Buican  2005,  2006).  Because  of  the  large 
number  of  optical  interfaces  in  the  PEM  stack,  the  light 
throughput  of  these  instruments  is  relatively  low. 

Third  generation  PEM/FT  spectrometers  (Fig.  lc) 
eliminate  the  high  light  loss  associated  with  the  many 
optical  interfaces  in  the  stack  by  replacing  the  latter  with  a 
“virtual”  stack  (VS):  this  is  a  solid  bar  of  optical  material 


that  has  only  two  optical  interfaces  but  which,  in  terms  of 
elastic  wave  propagation,  behaves  as  if  divided  into 
stacked,  isolated  regions.  The  VS  effect  is  created  by  a 
phased  array  (PA)  of  PZT  drivers  that  are  individually 
driven  at  well-determined  relative  amplitudes  and  phases. 
The  resulting  VSPA  concept  (Buican,  2008)  is  described 
below. 

3.  THE  VSPA  PHOTOELASTIC  MODULATOR 
3.1  Technical  Context 

It  was  shown  by  Canit  and  Badoz  (Canit  and  Badoz, 
1983)  that  small,  low-power  piezoelectric  transducers 
attached  to  the  sides  of  a  bar-shaped  elastic  solid  excite 
stress-waves  propagating  along  the  long  axis  of  the  bar 
medium  (Fig.  2).  Canit  and  Badoz  demonstrated  stress 
amplitudes  adequate  to  generate  half-wave  retardation  in 
light  propagating  perpendicular  to  the  bar  by  inducing  a 
standing  elastic  wave  along  the  bar’s  length,  rather  than 
across  the  medium  as  is  the  case  in  Kemp-type  photoelas¬ 
tic  modulators  (Kemp,  1969).  The  Canit  and  Badoz 
(C&B) 


Fig.  2.  Two  important  PEM  designs.  In  the  Kemp  design  (top),  a  funda¬ 
mental  mode  standing  elastic  wave  is  set  up  perpendicular  to  the  direc¬ 
tion  of  the  polarization-modulated  light  beam,  and  the  PZT  oscillations 
are  collinear  with  the  direction  of  the  standing  wave.  (Standing  elastic 
waves  are  shown  here  by  ribbon-like  representations  of  their  amplitude 
envelopes;  white  double-ended  arrows  represent  the  direction  of  propa¬ 
gation  of  the  elastic  waves;  and  light  beams  are  represented  by  green 
arrows.)  The  C&B  design  (bottom)  is  also  based  on  standing  elastic 
waves  perpendicular  to  the  direction  of  the  light  beam,  but  with  PZTs 
mounted  perpendicular  to  the  standing  wave.  For  longer  bars  (bottom 
right),  the  driving  frequency  is  tuned  to  an  overtone  of  the  axial  reson¬ 
ance.  In  this  case,  adjacent  maxima  (antinodes)  of  the  standing  wave 
have  opposite  phase  (represented  here  by  the  red  and  purple  colors  of  the 
amplitude  envelopes). 

design  dispenses  with  the  need  for  PZTs  and  optical  ele¬ 
ments  to  be  finely  tuned  to  a  common  resonant  frequency; 
this  simplifies  the  manufacturing  process  and  reduces 
manufacturing  cost.  The  ceramic  non-resonant  PZTs  used 
in  the  C&B  design  typically  require  the  use  of  a  hard  ad¬ 
hesive  for  good  mechanical  coupling,  which  can  give  rise 
to  static  stresses  and  thus  to  parasitic  static  birefringence. 
By  exciting  overtones  in  longer  bars  (Fig.2,  bottom  right), 
Canit  and  Badoz  showed  that  the  light  beam  can  be 
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moved  away  from  the  PZTs  and  the  associated  static  bire¬ 
fringence. 

The  C&B  PEM  design  is  interesting  because  its  re¬ 
tardation  amplitude  scales  with  bar  length  and  number  of 
attached  PZTs.  This  relationship  is  the  consequence  of  the 
PZTs  being  mounted  along  the  standing  wave  rather  than 
at  its  ends,  as  in  the  Kemp  design.  Consequently,  mechan¬ 
ical  power  can  be  injected  into  the  optical  element  at  each 
antinode  and,  except  for  losses  due  to  attenuation,  it  even¬ 
tually  reaches  the  light  beam  and  contributes  to  the  local 
birefringence  amplitude.  The  limitation  of  the  design 
comes  from  the  fact  that  the  length  of  the  optical  path — 
one  of  the  transverse  dimensions  of  the  bar — does  not 
change  as  the  length  of  the  bar  and  number  of  PZTs  in¬ 
crease. 

The  VSPA  concept  preserves  the  use  of  an  elastic 
standing  wave  along  the  length  of  a  bar-shaped  optical 
element  and  of  an  array  of  PZTs  placed  along  that  wave. 
However,  instead  of  using  the  end  surfaces  of  the  bar  to 
define  mechanical  resonance,  we  use  active  control  of  the 
individual  PZT  amplitude  and  phase  in  order  to  isolate 
adjacent  antinodes  and  drive  them  all  in  phase.  This  al¬ 
lows  us  to  have  an  optical  path  that  is  coaxial  with  the 
elastic  standing  waves.  Consequently,  in  a  VSPA  device, 
the  optical  path  length  and  total  retardation  amplitude 
scale  with  bar  length  while  the  local  stress  amplitude  does 
not. 

3.2  Phased  Arrays  of  PZTs 

We  introduce  the  conceptual  framework  behind  the  VSPA 
technology  by  means  of  a  simple  scalar  1-D  model  of 
elastic  wave  excitation  and  propagation  in  a  bar-shaped 
optical  element  with  attached  PZTs.  We  represent  the 
latter  by  point  sources.  The  simple  model  also  assumes 
that  the  stress  wave  propagates  along  one  of  the  principal 
axes  of  the  stress  tensor;  that  the  other  two  axes  do  not 
change  direction  over  the  domain  of  propagation  of  the 
elastic  waves;  and  that  the  propagation  medium  at  equili¬ 
brium  is  isotropic.  Consequently,  the  principal  axes  of 
stress,  stress-induced  dielectric,  and  refractive  index,  ten¬ 
sors  coincide  (Born  and  Wolf,  1970)  and  both  stress  and 
induced  birefringence  can  be  treated  as  scalars  that  are 
proportional  to  each  other.  Fig.  3a  shows  a  point  source 
L(X't)  A  0 

- • - x 


L(x,t)  R(x,t) 


Fig.  3.  (a)  A  point  source  at  x  =  x0  generates  a  local  oscillation  of  ampli¬ 
tude  A  o  and  phase  cp0  which  propagates  to  the  left  and  right  as  waves 
L(x,t)  and  R(x,t),  respectively,  (b)  A  set  of  N  point  sources  contained 
within  a  finite  interval  and  the  waves  L(x,t)  and  R(x,t)  propagating  out¬ 
side  the  interval. 

at  x0  which  generates  a  stress  (and  associated  birefrin¬ 
gence)  wave  that  propagates  along  the  x  axis  and  gene¬ 


rates  a  local  oscillation  of  amplitude  A0  (x) ,  phase  (x) , 
and  angular  frequency  co.  If  the  attenuation  coefficient  in 
the  medium  is  /?,  waves  propagating  to  the  left  and  right 
of  the  source  are  given  by  L(x,t)  -  A0el{cot+Kix~Xo))  and 

R(x,t)  =  where  A0  =  A0el*°  is  the 

complex  amplitude  describing  the  oscillations  at  the 
source;  and  k  is  the  complex  wavenumber  k  =  k-ij3 , 
where  k  =  In  /  X  is  the  (real)  wavenumber. 

Fig.  3b  shows  a  set  of  N  point  sources  at 
Xj  <x2  <...<xN  and  oscillating  with  the  same  angular 
frequency  co  and  generating  N  sets  of  waves  having  am¬ 
plitudes  and  phases^. and^.  at  their  respective  sources. 

The  sources  are  contained  within  a  finite  interval  along 
the  x  axis,  with  L(x,t)  and  R(x,  t)  being,  respectively,  the 
waves  propagating  to  the  left  and  to  the  right  of  that  inter¬ 
val.  The  left-  and  right-propagating  waves  are  given  by 


L(x,t)  = 


R(x,t )  = 


ZAA'K 

V  j=]  J 


and 


(la) 

(lb) 


The  complex  amplitudes  of  the  two  waves  are  seen  to  be 
essentially  the  Fourier  transforms  of  the  spatial  distribu¬ 
tion  of  sources,  evaluated  at  ±k.  Thus,  the  amplitudes  of 
the  two  waves  are  constant  outside  the  source  interval, 
except  for  the  effect  of  attenuation  in  the  medium.  Conse¬ 
quently,  if  the  amplitude  of  one  of  these  waves  vanishes 
on  the  boundary  of  the  source  interval,  than  its  amplitude 
is  zero  anywhere  outside  the  boundary  and  the  purely 
geometric  boundary  point  where  the  amplitude  vanishes 
behaves  like  a  mirror  for  the  waves  generated  by  the  set 
of  sources. 


3.3  Active  Reflectors 

The  waves  propagating  out  of  the  source  interval  va¬ 
nish  if  the  wave  complex  amplitudes  in  Equations  (la) 

N  _  N 

and  (lb)  vanish,  ^A je~lKXj  =0  and  ^A jelKXj  =0. 

7=1  7=1 

The  above  equations  are  thus  the  wave  cancellation  condi¬ 
tions.  If  only  one  of  these  conditions  is  satisfied,  then 
there  is  no  wave  propagation  on  the  corresponding  side  of 
the  source  interval.  The  above  equations  each  have  non¬ 
trivial  solutions  for  A>  2;  in  the  simplest  case,  N=  2,  one 
source  can  always  be  driven  so  as  to  cancel  the  wave  ar¬ 
riving  from  the  other.  In  this  one-sided  wave  confinement, 
the  cancelling  source  behaves  as  a  virtual,  active,  reflec¬ 
tor.  The  complex  amplitudes  for  the  virtual  reflectors  at 
Xj  and  x7V(Fig.  3b)  can  easily  be  obtained  by  solving  the 

wave  cancellation  equations  for  Aj  and  respective¬ 
ly.  These  solutions,  which  depend  in  general  on  the  posi- 
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tions  of  all  sources,  directly  provide  a  global,  open-loop 
approach  to  driving  an  active  reflector.  This  can  be  quite 
onerous  in  terms  of  sensors  and  data  processing  require¬ 
ments  and  thus  impractical.  However,  the  existence  of  a 
nontrivial  solution  to  each  complex  amplitude  equation 
makes  possible  a  local  approach  to  driving  an  active  ref¬ 
lector,  in  which  a  local  sensor  or  small  set  of  sensors,  to¬ 
gether  with  a  closed-loop  controller,  keep  the  wave  ampli¬ 
tude  at  the  active  reflector  at  a  certain  acceptable  error 
level.  Interestingly,  this  approach  to  noise  cancellation 
dates  back  to  the  1930s  (Lueg,  1934;  see  also  Tokhi  and 
Veres,  2002). 

An  active  reflector  can  also  be  used  as  an  isolator 
that  separates  two  1-D  domains  and  prevents  interactions 
between  them  through  the  cross-propagation  of  waves.  As 
before,  consider  a  set  of  M  point  sources  to  the  left  of  the 
source  to  be  used  as  isolator,  positioned  at  xtj  such  that 

xn  <  xn  <  •••  <  xim  (Fig-  4)  and  generating  M  sets  of  waves 
L(x,t)  ;  R(x,t) 


Fig.  4.  If  driven  with  the  right  phase  and  amplitude  (see  text),  the  point 
source  at  x0  can  block  the  propagation  of  waves  from  one  side  to  the 
other  and  thus  effectively  isolate  the  two  sets  of  sources. 

with  amplitudes  and  phases  Atj  and  (/>tj  at  their  respective 

sources.  Similarly,  consider  a  set  of  N  point  sources  to  the 
right  of  the  isolator,  and  use  corresponding  parameters 
having  the  index  “r”  instead  of  “1”.  From  the  wave  cancel¬ 
lation  equations,  the  waves  generated  by  the  sources  to 
the  right  of  the  isolator  do  not  propagate  through  it 
(L(x,t)= 0  in  Fig.  4)  if  the  isolator  source  is  driven  with 
the  complex  amplitude  Ar  given  by 

N 

Ar  =  Arje~lK{Xi~Xo)) ,  where  x0is  the  position  of  the 

M 

isolator  source.  Similarly,  the  complex  amplitude  for 
which  the  isolator  blocks  the  propagation  of  the  waves 
generated  by  the  sources  to  its  left  (R(x,t)= 0  in  Fig.  4), 

M 

At  -  A/y-gIK(**~*o)) .  Under  the  assumption  of  linearity, 

7=1 

driving  the  isolator  with  a  complex  amplitude  equal  to  the 
sum  A0  =  Ar  +  At  will  isolate  the  two  sets  of  sources. 

(See  Fig.  5;  also  note  the  simple  example  of  active  isola¬ 
tor  shown  in  Fig.  A2c  of  the  Appendix.) 

Fig.  5  shows  how  adjacent  standing  wave  antinodes 
can  be  decoupled  by  means  of  an  active  isolator.  In  a 
thought  experiment,  two  C&B  fundamental-frequency 
PEMs  are  physically  stacked  (top  right)  and  then  brought 
together  until  they  fuse  (center  right)  and  form  one  single 
bar  of  twice  the  original  resonator  length.  This  is  similar 
to  a  first-overtone  C&B  PEM  (Canit  and  Badoz,  1983),  in 
which  the  interfering  elastic  waves  propagate  through  the 
central  node  and  thus  cause  the  two  antinodes  to  oscillate 
with  opposite  phase.  By  adding  a  node  driver  (bottom 


right)  oscillating  with  the  appropriate  relative  amplitude 
and  phase,  the  two  antinodes  become  independent  again 
and  can  be  driven  in  phase  just  like  the  separate  PEMs  in 
the  physical  stack  (top  right).  Thus,  the  bar-shaped  device 
in  the  bottom-right  diagram  can  be  regarded  as  a  “virtual” 


Fig.5.  Conceptual  path  from  a  stack  of  PEMs  to  a  Virtual-Stack  PEM. 
Fundamental-mode  C&B  PEMs  (left)  can  be  physically  stacked  (top 
right)  to  increase  total  retardation  amplitude,  but  at  the  cost  of  reduced 
light  throughput  because  of  reflection  losses  at  the  multiple  optical  inter¬ 
faces.  If  the  PEMs  are  brought  into  contact  with  each  other  so  as  to  fuse 
into  a  single  bar-shaped  PEM  (center  right),  additional  light  losses  are 
eliminated  but  phase  reversal  at  the  central  node  causes  the  retardations 
in  adjacent  antinodes  to  cancel  each  other  out.  The  addition  of  a  “node 
driver”  (bottom  right)  and  active  vibration  cancellation  isolates  the  adja¬ 
cent  antinodes  which  can  now  oscillate  in  phase.  (See  Fig.  A2c  for  the 
numerical  model  of  a  similar  configuration.) 

stack  driven  by  the  phased  array  consisting  of  one  node 
and  two  antinode  drivers — the  simplest  example  of  a 
VSPA  PEM. 


3.4  Active  Resonance  Cells 

If  both  wave  amplitudes  in  Equations  (la)  and  (lb)  va¬ 
nish,  there  is  active  confinement  of  the  waves  within  the 
source  interval.  A  nontrivial  solution  always  exists  for  N 
>  3.  In  the  simplest  case,  N=  3  (Fig.  6,  center),  the 

A 

Pt 


Nan 


Fig.  6.  Alternative  conceptual  path  from  a  C&B  PEM  to  a  Virtual-Stack 
PEM.  The  equivalent  of  a  fundamental-mode  C&B  PEM  (top)  can  be 
implemented  by  means  of  an  ARC  in  a  bar-shaped  optical  element  (cen¬ 
ter)  (see  also  Fig.  A2).  Adjacent  ARCs  can  share  node  drivers  and  a 
virtual  stack  of  C&B -like  PEMs  can  be  implemented  within  a  bar  (see 
also  Fig.  A2c).  (A — antinode  driver;  N — node  driver.  PZT  drivers  di¬ 
rectly  above  each  other  are  of  the  same  type.) 
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two  outside  sources  can  be  seen  as  actively  confining  a 
standing  wave  inside  a  virtual  resonator.  The  standing 
wave  appears  to  be  driven  by  the  middle  source.  For  this 
reason,  we  refer  to  the  two  outside  sources  in  this  confi¬ 
guration  as  “node  drivers” — they  are  placed  at  nodal 
points  of  the  confined  standing  wave — and  to  the  middle 
source  as  an  “antinode  driver”,  even  though  it  need  not  be 
positioned  exactly  at  a  maximum  of  the  wave  amplitude. 
In  reality,  all  three  sources  drive  the  confined  standing 
wave. 

For  N=3,  the  optical  medium  and  PZT  driver  array 
shown  in  Fig.  6  (center)  form  a  “virtual”  resonator  that, 
without  being  bounded  by  physical  discontinuities,  be¬ 
haves  like  the  physical  resonator  shown  in  Fig.  6  (top)  in 
terms  of  wave  confinement.  Satisfying  both  equations  for 
the  complex  amplitudes  means  keeping  the  driving  ampli¬ 
tudes  and  phases  of  the  PZTs  in  certain  relationships — the 
PZTs,  together  with  their  drive  and  control  electronics, 
thus  form  a  phased  array.  This  phased  array  and  the  me¬ 
dium  in  which  it  excites  elastic  waves  form  what  we  call 
an  Active  Resonance  Cell  (ARC). 

The  lack  of  physical  discontinuities  in  the  ARC 
means  that  the  propagation  of  optical  waves  along  the 
optical  axis  of  the  device  is  not  affected  by  scattering  and 
absorption  losses  normally  caused  by  such  discontinuities. 
This  is  a  fundamental  advantage  if  achieving  high  retarda¬ 
tion  amplitudes  requires  the  use  of  multiple  resonators. 

The  effectiveness  of  an  ARC  is  described  by  its  ef¬ 
fective  stress  amplitude  (see  the  Appendix).  The  value  of 
the  effective  amplitude  is  given  for  each  of  the  examples 
shown  in  the  Appendix.  These  values  illustrate  the  fact 
that  local  stresses  and  birefringences  add  up  constructive¬ 
ly  along  an  ARC’S  optical  path  and  yield  a  high  retarda¬ 
tion  amplitude.  This  is  due  to  the  essentially  standing  na¬ 
ture  of  the  waves  inside  the  ARCs,  as  well  as  to  the  ab¬ 
sence  of  phase  reversal  at  node  drivers  (isolators). 

ARCs  are  robust.  As  the  examples  in  the  Appendix 
show  (Fig.  A2),  wave  confinement  can  be  established 
even  when  the  geometric  and  elastic  wave  propagation 
parameters  are  far  removed  from  those  required  in  a  phys¬ 
ical  resonator.  The  effective  stress  amplitudes  also  re¬ 
main  high  under  such  adverse  conditions.  Thus,  both 
manufacturing  variability  and  thermal  effects  can  be  easi¬ 
ly  compensated  for  electronically  during  operation  by 
adjusting  PZT  amplitudes  and  phases. 

4.  MODELING  AND  COMPUTATION 
4.1  Finite  Element  Modeling 

The  Finite  Element  Method  (FEM)  (see  Pepper  and 
Heinrich,  2006,  for  a  good  introduction  to  FEM)  is  a  nu¬ 
merical  technique  for  solving  partial  differential  equations 
by  discretizing  their  spatial  domain  into  a  mesh  of  small 
but  arbitrarily  shaped  regions  (finite  elements).  This  re¬ 
sults  in  matrix  equations  that  relate  the  input  at  certain 


points  (nodes)  in  the  elements  to  outputs  at  the  same 
points.  A  short  history  of  the  method  can  be  found  in 
Pepper  and  Heinrich  (Pepper  and  Heinrich,  2006).  Several 
well-established  commercial  packages  based  on  FEM 
algorithms  are  available.  For  this  work,  we  chose  COM- 
SOL  Multiphysics  (COMSOL,  2006)  because  of  its  ease 
of  use  and  flexibility  in  incorporating  and  integrating 
models  for  a  multitude  of  interacting  physical  phenomena. 

The  first  published  structural  mechanics  computer  model¬ 
ing  of  PEMs  that  we  are  aware  of  came  from  Canit  and 
coworkers  (Yang  et  al.,  1995),  who  used  a  commercial 
FEM  package  to  study  the  longitudinal  elastic  eigen- 
modes  of  their  bar-shaped  PEMs.  Figures  7a  and  7b  show 
plots  of  the  first  two  eigenmodes  corresponding  to  C&M 
PEM  operation.  These  plots  were  generated  as  part  of  our 
own  survey  of  bar  eigenmodes. 


Fig.  7.  COMSOL  models  of  the  bar-shaped  PEM.  First  (a)  and  second 
(b)  eigenmodes  of  a  rectangular  bar  (see  the  text  for  details);  An  impor¬ 
tant  higher  eigenmode  of  the  bar:  boundary  deformation  and  color-coded 
vertical  displacement  (c)  and  principal  strain  vectors  (d);  Virtual  reflec¬ 
tor  effect  showing  a  vertically  polarized  shear  wave  preferentially  prop¬ 
agating  leftwards  (A — antinode  driver  PZTs;  N — node  driver  PZTs). 
The  eigenmode  illustrations  show  only  a  quarter  of  the  bar  (lower  right 
quadrant  of  the  cross-section)  which  is  assumed  to  be  symmetrical. 

The  C&B-type  eigenmodes  are  dominated  by  longi¬ 
tudinal  (compression)  waves  that  are  reflected  at  the  bar’s 
ends.  As  can  be  gleaned  from  Fig.  7a  and  Fig.  7b,  these 
are  accompanied  by  transverse  compression  that  is  the 
same  in  both  transverse  directions  and  therefore  does  not 
result  in  significant  birefringence  in  a  transverse  cross- 
section  of  the  bar.  Consequently,  light  traveling  along  the 
length  of  the  bar  as  in  a  VSPA  PEM  does  not  record  sig¬ 
nificant  retardation  between  the  two  directions  of  polari¬ 
zation,  while  light  traveling  as  in  a  C&B  PEM  does. 

At  the  other  extreme,  a  longitudinally  nonpropagating 
mode  (Fig.  7c)  similar  to  that  used  in  octagonal  Kemp- 
style  PEMs  (Hinds  Inc.,  Hillsboro,  Oregon)  can  be  estab¬ 
lished  (Fig.  2,  top)  by  reflection  of  compression  waves 
propagating  perpendicular  to  the  length  of  the  bar.  As  in 
the  conventional  Kemp  design,  this  requires  precise  grind¬ 
ing  of  the  side  surfaces,  which  furthermore  must  apply  to 
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the  full  length  of  the  bar.  This  is,  however,  a  highly  effec¬ 
tive  way  of  generating  high  birefringence  along  the  length 
of  the  bar  since  the  stress  and  refractive  index  tensors 
have  the  eigenvectors  for  the  highest  and  lowest  eigenva¬ 
lues  in  the  transverse  cross-section  of  the  bar  (see  Fig.  7  d 
for  the  corresponding  strain  eigenvectors).  Furthermore, 
this  high  birefringence  can  extend  to  most  of  the  length  of 
the  bar. 

We  also  use  COMSOL  Multiphysics  to  generate 
time-dependent  simulations  showing  the  propagation  of 
PZT-generated  elastic  waves.  For  example,  Fig.  7e  shows 
an  active  mirror  (node  driver)  as  described  in  Section  3 
that  reflects  (with  some  losses)  a  pure  shear  waves  propa¬ 
gating  from  a  source  (antinode  driver)  to  the  left  of  the 
reflector.  This  is  important  because  the  compromise  be¬ 
tween  extreme  manufacturing  precision  and  low  PEM 
performance  can  be  found  among  the  many  propagating 
modes  characteristic  of  the  solid  bar  (Auld,  1990)  in  con¬ 
junction  with  VSPA  driving  and  control  technology,  and 
these  modes  generally  exhibit  significant  shear. 

4.2  Computation 

We  are  currently  using  both  free  open-source  and 
commercial  off-the-shelf  (COTS)  finite  element  modeling 
(FEM)  software  to  create  comprehensive  multiphysics 
models  of  the  VSPA  PEM  optical  element.  The  models 
are  developed  and  tested  on  HERD-0,  a 
small  in-house  PC  cluster  (Fig.  8). 

The  HERD-0  cluster  uses  many 
open  source  software  packages  such  as 
the  Rocks  cluster  management  system 
(Papadopoulos  et  al.,  2001)  and  the 
Trilinos  library  of  object-oriented  linear 

Fig.  8.  The  HERD-0  PC  cluster  assembled  at  SEA 
LLC  for  physical  modeling  work.  Each  of  the 
nine  nodes  has  a  dual  core  Athlon  64  processor 
and  3  GB  of  RAM.  The  nine  PCs  are  connected 
by  a  private  Gigabit  Ethernet  network  and  are 
powered  by  uninterruptible  power  supplies. 
HERD-0  exceeded  56  Gflops/s. 

algebra  tools  and  solvers  (Heroux  et  al.,  2003).  The  com¬ 
mercial  COMSOL  Multiphysics  package  (COMSOL  AB, 
2006)  is  installed  on  the  cluster’s  head  node  and  is  used 
for  multiphysics  modeling  of  the  interacting  mechanical, 
thermal  and  optical  processes  in  the  PEM  bar  and  PZT 
array. 

An  interface  between  COMSOL  and  the  Trilinos 
solvers  running  in  parallel  on  the  cluster  was  developed  in 
house;  this  interface  allows  the  problem  set  up  in  COM¬ 
SOL  to  be  transferred,  via  COMSOL  Script  (COMSOL 
AB,  2006),  Java/JNI  and  C++  code,  to  the  appropriate 
cluster-based  Trilinos  solvers,  and  the  results  transferred 
back  to  COMSOL  for  postprocessing  and  display.  A  large 
PC  cluster  at  the  Army  Research  Laboratory  (ARL)  Major 
Shared  Resource  Center  (MSRC),  Aberdeen  Proving 
Ground,  MD,  will  be  used  for  production  runs  once  the 


FEM  model  has  been  debugged  and  tested  on  the  local 
cluster. 

5.  PRELIMINARY  DESIGN  CONSIDERATIONS 

The  diagrams  in  Figures  1,  5  and  6  illustrate  the  use 
of  two  nonresonant  ceramic  PZTs  similar  to  the  ones  used 
by  Canit  et  al.  (Fig.  2  and  Canit  and  Badoz,  1983;  Yang  et 
al.,  1995)  for  each  node  and  antinode  driver.  These  PZTs 
are  placed  on  opposite  faces  of  the  bar-shaped  optical 
element.  Additional  configurations  consisting  of  four  or 
more  PZTs  per  node  or  antinode  driver  are  under  consid¬ 
eration  and  seem  to  be  required  for  exciting  vibration 
modes  with  high  associated  birefringence  (see  Figures  7c 
and  7d  for  an  example  of  such  a  mode).  The  driving  and 
control  electronics  for  the  phased  array  of  PZTs  deployed 
along  the  length  of  the  bar  may  need  to  be  supplemented 
by  additional  phase  controls  between  PZT  on  different 
faces  of  the  bar.  These  and  other  recent  work  on  the  effi¬ 
cient  driving  of  the  phased  array  will  be  described  else¬ 
where. 

The  use  of  relatively  inexpensive  nonresonant  ceram¬ 
ic  PZTs  will  result  in  a  considerable  reduction  in  the  cost 
of  high  retardation  amplitude  PEMs.  The  elimination  of 
the  numerous  optical  interfaces  present  in  a  PEM  stack 
will  also  considerably  reduce  manufacturing  costs.  Fur¬ 
thermore,  the  use  of  low-profile  PZTs  with  a  large  contact 
area  will  likely  result  in  rugged  designs  suitable  for  use  in 
the  field  and  on  mobile  platforms  such  as  combat  vehicles 
and  helicopters. 

6.  DEFENSE  APPLICATIONS 

We  have  previously  discussed  (Buican  and  Carrieri, 
2004;  Buican  and  Carrieri,  2006)  a  range  of  explosive 
trace  detection  and  CB  defense  applications  that  are  made 
possible  by,  or  can  benefit  considerably  from,  the  use  of 
PEM/FT  spectrometers.  These  applications  range  from 
novel  spectrometric  techniques  for  the  detection  and  rec¬ 
ognition  of  surface  contaminants  (Buican  and  Carrieri, 
2006;  Carrieri  et  al,  1999)  to  panoramic  hyperspectral  IR 
imagers  for  detecting  and  tracking  chemical  clouds  (Car¬ 
rieri,  2003). 

Any  application  that  requires  the  speed,  broadband 
characteristics  and  light  throughput  of  a  PEM/FT  spec¬ 
trometer  would  benefit  even  more  from  the  third- 
generation  VSPA  PEM  technology.  As  previously  dis¬ 
cussed,  the  latter  improves  on  the  light  throughput  of 
stacked-PEM  instruments  while  allowing  for  more  rugged 
and  compact  instrument  designs  that  are  better  suited  for 
field  use.  In  fact,  because  of  their  inherent  advantages 
over  older  approaches,  VSPA  PEMs  will  likely  be  the 
common  enabling  technical  element  for  these  novel  detec¬ 
tion  technologies.  As  these  techniques  are  strategically 
important  by  virtue  of  their  likely  impact  in  combat  zones 
and  in  homeland  security,  VSPA  PEMs  themselves  be¬ 
come  the  nexus  of  a  strategic  advantage. 
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We  briefly  review  in  what  follows  three  strategically 
important  spectral  techniques  for  use  in  explosive  trace 
detection  and  CB  defense  and  for  which  VSPA  will  likely 
turn  out  to  be  the  enabling  technology. 

We  described  multivariate  PEM/FT  spectrometry,  M- 
PEMS  (Buican  and  Carrieri,  2006),  as  an  active  spectral 
technique  that  provides  correlated  excitation,  emission 
and  lifetime  fluorescence  spectra.  These  three- 
dimensional  spectra  are  inherently  more  likely  to  provide 
the  basis  for  reliable  discrimination  than  regular  spectral 
data.  Furthermore,  multivariate  spectrometers  could  be 
used  in  standoff  configurations  to  analyze  small  sample 
volumes,  including  small  volumes  of  surface  contami¬ 
nants.  And,  when  used  in  conjunction  with  ultrashort 
pulse  lasers,  multivariate  spectrometers  could  analyze 
trace  surface  contaminants  of  low  vapor  pressure  and/or 
low  intrinsic  fluorescence. 

Thermal  luminescence  (TL)  is  a  technology  aimed  at 
solving  the  standoff  detection  problem  of  contaminated 
water  reservoirs  and  terrestrial  surfaces  (Carrieri  et  al., 
1999).  In  TL,  the  energy  of  an  irradiating  beam  is  ab¬ 
sorbed  into  the  target  material;  by  scanning  the  irradiated 
surface,  the  PEM/FTIR  spectrometer  produces  high-speed 
interferograms  that  are  transformed  into  graybody  spectra 
of  the  TL.  During  beam  irradiation,  a  thermal  gradient  is 
produced  at  the  target  surface;  inside  the  thermal  nonequi¬ 
librium  heating  period,  a  detection  window  of  opportunity 
develops  and  the  absorption  bands  of  contaminant  are 
revealed.  These  bands  are  quite  strong  provided  the  gray¬ 
body  envelope  shifts  at  maximum  rate  (i.e.,  the  maximum 
gradient  condition),  which  requires  high  interferogram 
acquisition  rates  that  can  only  be  provided  by  PEM/FTIR 
spectrometers. 

PANSPEC  (Panoramic  Infrared-Imaging  Spectro- 
radiometer)  is  a  computer-optimized  spectroradiometer 
design  for  monitoring  a  panoramic  infrared  environment 
for  the  presence  of  chemical  clouds,  and  for  tracking  such 
clouds  (Carrieri,  1997a,  1997b,  2003).  It  combines  imag¬ 
ing  PEM/FTIR  and  photopolarimeter  technologies  and 
integrates  a  wide-angle  camera  in  forming  a  unified  pas¬ 
sive  imaging  remote  sensor  and  active  laser  beacon  sys¬ 
tem.  PANSPEC  provide  a  standoff  detection  capability 
for  disseminated  chemical  weapons;  specifically,  (1)  de¬ 
tection  of  gaseous  chemical  warfare  agents  in  the  atmos¬ 
phere  to  a  range  of  approximately  50  km  (a  figure  of  me¬ 
rit  for  prototype  construction);  (2)  tracking  of  the  targeted 
hazardous  cloud  mass  over  a  panoramic  field  of  view;  and 
(3)  omnidirectional  communications  of  cloud  presence 
and  heading  by  means  of  an  encrypted  laser  beam  carrier, 
projected  into  the  exact  panorama  from  which  imaging  is 
done. 


CONCLUSIONS 

The  VSPA  concept  shows  considerable  promise  as 
the  basis  for  high  retardation  amplitude  PEMs  for  use  in 


ultra-high  performance  FT  spectrometers.  While  addition¬ 
al  modeling  and  experimental  work  is  required  in  order  to 
refine  the  concept  and  demonstrate  it  experimentally,  we 
believe  that  VSPA  PEM/FT  spectrometry  has  the  poten¬ 
tial  to  become  the  enabling  technology  for  fast,  sensitive, 
and  reliable  optical  detection  and  identification  of  specific 
chemical  and  biological  agents  and  contaminants. 

APPENDIX 


We  define  the  effective  amplitude  and  phase  for  the 
parameters  of  interest — in  this  case,  stress  and  stress- 
induced  birefringence — as  the  (real)  amplitude  and  phase 
of  the  average  value  of  the  complex  amplitude  over  the 
optical  path  through  the  ARC  or  well-defined  portion 


thereof,  A  -  A  ffel(Peff 

’  eff  eff 


A(x,t)dx,  where  A 


is  the 


interval  of  interest  inside  the  ARC  and  A(x,  t)  is  the  com¬ 
plex  amplitude  of  the  oscillations  w(x,t)  inside  the  inter¬ 
val,  w(x,t)  =  A{x,t)elcot . 


The  examples  given  below  show  the  amplitude  and 
phase  distributions  for  ARCs  with  various  geometries  and 
elastic  wave  propagation  properties.  The  antinode  driver 
is  positioned  at  the  origin.  The  wave’s  peak  amplitude  is 
set  to  one  amplitude  unit  and  its  phase  at  the  origin  to 
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Fig.  Al.  The  three  sources  forming  an  ARC  (top)  in  the  examples  shown 
in  Fig.  A2  and  the  equivalent  conventional  resonator  bounded  by  two 
reflecting  surfaces  (bottom).  Of  course,  the  conventional  resonator  ade¬ 
quately  functions  only  for  reflector  separations  that  are  integer  multiples 
of  X/2.  The  node  drivers  are  situated  at  jti,  jc3,  and  the  antinode  driver  is 
at  the  origin. 


zero.  Please  note  that  the  waves  in  all  the  examples  below 
are  fully  contained  within  the  ARC  despite  large  devia¬ 
tions  from  standing-wave  conditions. 


=  0.25;  p  =  0;  Aeff=  0.636 


(b)  X!  = -0.4,  x2=0,x3=  0.2;  p  =  0.5;  Aeff=  0.596 
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(c)  xi  =  -0.5,  x2  =  0,  x3  =  0.5;  (3  =  0;  Aeff  =  0.636 
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(d)  xi  = -0.1,  x2=0,x3=  0.2;  p  =  0.5;  Aeff  =0.557 


Fig.  A2.  Each  row  of  drawings  shows  the  amplitude  (left)  and  phase 
(right)  profile  for  an  ARC  (black  lines).  The  red  vertical  bars  show  re¬ 
spectively  the  values  of  the  amplitudes  and  phases  for  the  three  sources 
of  the  ARC.  The  dashed  lines  show  respectively  the  effective  amplitude 
and  phase  for  each  ARC  interval,  and  the  point-dashed  line  shows  the 
values  for  the  entire  ARC.  The  data  were  normalized  for  unit  peak  wave 
amplitude  and  for  zero  wave  phase  at  the  origin.  For  each  example  we 
give  the  positions  of  the  three  sources  (in  wavelengths);  the  attenuation 
coefficient  (3  (per  wavelength);  and  the  effective  amplitude  over  the 
entire  ARC.  (a)  The  node-antinode  separation  is  A/4  and  there  is  no 
attenuation;  this  corresponds  to  the  fundamental  resonance  and  is  a 
special  case  in  that  the  antinode  source  amplitude  is  zero  and  the  ARC  is 
reduced  to  only  two  sources.  The  effective  amplitude  is  0.636,  the  cha¬ 
racteristic  value  for  a  conventional  standing  wave  between  two  adjacent 
nodes.  The  node  drivers  are  characteristically  ahead  of  the  wave  phase  at 
the  origin  by  7i/2.  (b)  A  highly  nonsymmetric  ARC  with  high  attenua¬ 
tion;  although  there  is  some  change  in  wave  phase  within  the  ARC,  the 
effective  amplitude  is  still  high,  (c)  The  node-antinode  separation  is  A/2 
and  there  is  no  attenuation;  each  half  of  the  ARC  becomes  itself  an  ARC 
as  in  (a)  and  the  antinode  driver  can  be  regarded  as  an  isolator,  (d)  A 
highly  nonsymmetric  ARC  that  is  very  far  from  standing-wave  geome¬ 
try.  Adversity  notwithstanding,  the  effective  amplitude  is  still  87.6%  of 
the  value  for  the  ideal  standing  wave  in  (a). 
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